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ABSTRACT 

In this article, models are constructed of the compact jet in GRS 1915+105 dur- 
ing an epoch of optimal data capture. On April 02, 2003, the object was observed 
in the hard X-ray/soft gamma ray band (INTEGRAL), hard X-ray band (RXTE), 
near IR (ESO/New Technology Telescope) and the VLBA (8.3 GHz and 15 GHz). 
The source was in a so-called "high plateau state." The large radio flux provides 
high signal to noise ratios in the radio images. Thus, one can image the jet out to 
large distances (> 10^^ cm). This combined with the broadband coverage make this 
epoch the best suited f or modeling the jet. The paran i etric method developed in th e 
papers iGhisellini et all (|1985[ ): iGhisellini and Maraschil (fT989t ) : iGhisellini et all (|1996[ ): 
ISambruna et all (|1997D that has been successfully utilized in the realm of extragalactic 
radio jets is implemented. The basic model is one where external inverse Compton 
(EIC) scattering of accretion disk photons by jet plasma provides the hard X-ray 
powerlaw. Unlike AGN jets, it is found that the radio jet must be highly stratified 
in the transverse direction in order to produce the observed surface brightness dis- 
tribution in the radio images. Various jet models are considered. The jet power is 
Q « 3 — 4 X 10'^® ergs/sec if the hard X-ray powerlaw luminosity is from EIC in the 
jet and Q « 2 — 9 x 10'^'' ergs/sec if the X-rays are emitted from the accretion disk 
corona. These estimates indicate that the jet power can be as high as 60% of the total 
X-ray luminosity. 
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Galactic black hole accretion systems might serve as 
laboratories for AGN (active galactic nuclei) on highly 
contracted time scales. The black hole candidate, GRS 
1915-1-105, routinely ejects powerful radio emitting plasma 
out to large distances from the central black hole at rel- 
ativistic speeds and might be the Galactic analog of ra- 
dio l o ud AGN jRodriguez and Mirabell 1 19991 : iFender et all 
1 19991 : iRushton et aT |2010a|). These major flares are rela- 
tively rare compared to the so-called compact jet state 
l|Rushton et a l 2010a). A compact, optically thick, elongated 
structure i s frequently observed with long baseline inter- 
ferom etry jPhawan et alll200(]| : iFuchs et alll2003l : iRibo et all 
bOO^ V The low frequency emission is significantly less than 
in a major fiare. It is not clear if there is an analogous struc- 
ture in the context of AGN. The total power of the jet is also 
unclear. In this article, the analogy with AGN is pursued 
by modeling the jet u sing the standard AGN t echniq ue of 
IGhisellini etal (119851 ): IGhisellini and Maraschi' (19891) with 
the addition of external Compton s cattering (EIC) from 
(IGhisellini et "ai ll996l : ISambruna erai|[l997) ) . There is a cor- 
relation or more precisely a direct association of the for- 



mation of a compact jet with an excess of hard X-ray 
emission in the form of a power law (|Rushton et alllioiObI : 
iRodriguez et alll2008l ). The primary assumption of the mod- 
els of the jet presented below is t hat this X-ray emis sion 
is produced in the base of the jet (jMarkoff et alll2005l ). A 
scenario in which the X-rays are from Compton scattering 
in the accretion disk corona is also explored briefiy in this 
article with the same techniques. 

Since most of the luminosity is in the X-ray band, com- 
plete high energy coverage is mandatory for jet modeling. 
Furthermore, one needs the near IR photometry in multi- 
ple bands in order to constrain the synchrotron spectrum. 
Also, to break the degeneracy in the available models, one 
requires a faithful representation of the surface brightness 
distribution of the jet. It is a much simpler task to fit the 
jet spectrum if the distribution of emissivity is ignored. To 
meet this challenge in the presence of intrinsic source vari- 
ability, one needs simultaneous full coverage of the X-ray 
emission to high frequencies, IR photometry and sensitive 
VLBA imaging. The only simultaneous set of observations 
capable of fulfilling these requirements was the campaign 
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THE BASE OF THE NESTED SHELL JET MODEL 
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Figure 1. The shell model of transverse jet structure. 



presented in (|Fuchs et all 12001 12004| ). The strongest com- 
pact jet ever imaged witii VLBA was measured simultane- 
ously in tiie near IR (ESO/New Technology Telescope) and 
throughout the hard X-ray and soft gamma ray band (IN- 
TE GRAL and RXTE). For more details on the observations 
see iFuchs et all l|2003l. |2004| ) . 



considered as suggested in iDhawan et a 1 l|2000l ). Assuming 
P = 0.98 and 6 — 66° as was deduced from the kine- 
matics of the discrete ejections in iFender et all (|l999l ). it is 
found that the flux of blackbody radiation from the accre- 
tion disk is highly redshifted in the frame of reference of 
the jet plasma, so synchrotron self Compton emission (SSC) 
dominates the external Compton emission (EIC) and SSC 
is the source of the power law X-rays. Many models were 
constructed, but because the X-ray luminosity (of the un- 
resolved base of the bipolar jet) observed at earth is red- 
shifted by a factor of 5'*, the intri nsic X-ray luminosit y 
Lint(X-ray) ^ 3.5 x 10"° erg/sec (jLightman et allll975l '). 
Therefore the jet power, Q > 3.5 x 10*" erg/sec. This seems 
unreasonably large. Consequently, relativistic bulk motion 
is not considered in the following. These calculations indi- 
cate that a reasonable model should have at most mildly 
relativistic velocities that are consistent with the observed 
jet asymmetry (i.e., /3 « 0.4). EIC is therefore the source of 
the power law X-ray hard tail (not SSC) since the redshifted 
thermal X-ray flux is much larger than the local synchrotron 
flux within the base of the jet. This result is considered in 
the context of other estimates of the jet velocity in the lit- 
erature in the Discussion section. 



1 SOME PRELIMINARIES ON JET SPEED 

There is uncertainty in the bulk velocity of the compact 
jet. The Doppler enhancement of the intrinsic radiation is a 
critical parameter in the determination of the jet power, Q. 
The Doppler factor, 5, is given in terms of F, the Lorentz 
factor of the outflow; /?, the three velocity of the outflow 
and the angle of propagati on to the line o f sight, 9; 5 = 
l/[r(l - /3 cose)] (Lind an d Blandfordll 19851 ). For a jet, the 
apparent flux density, Sv is related to the intrinsic value by 
S^.obs = 5^+"5'^, int, where a is the spectral index (the con- 
vention adopted in this paper is SnU oc nu'"^'^^""). Assum- 
ing an intrinsically symmetric bipolar jet, the elevated flux 
of the approaching jet has been used to estimate S. Using a 
value of 66° < 6 < 70° that was de duced from discrete rel- 
ativist i c ejecta at ot h er ep ochs (see lMirabel and Rodriguej 
(| 19941 ): iFender et a] l| 19991 )). estimates of 0.1 < ^ < 0.54 
for the compact jet at various freque ncies and epochs of 
VLBI observatio ns have been obtained (iDhawan et a ll200d: 
iRibo et all [20041 ). For the epoch April 02, 2003, iRibo et all 
l|2004l) found disagreement between the jet velocity inferred 
from the flux asymmetries of the 15 GHz and 8.3 GHz VLBA 
images. They reconciled this by looking at the Clean Com- 
ponents, concluding that a value of P — 0.38 ± 0.04 was 
compatible with the data. Thus, /3 — 0.4 is the baseline 
value in the following calculations. 

The discussion above is based on direct application 
of the data to relativistic b eaming as was done for dis- 
crete ejecta in m ajor flares |Mirabel and Rodriguej 1 19941 : 
IFender et all 19991 ). In IDhawan et all (|20od ). the difference 
in the directly inferred speeds of the compact jet and the 
major flares was disconcerting. Thus, they conjectured on 
other explanations including undetec ted relativistic motion. 
Thus, before the iRibo etal l|2004l ) estimated value was 
adopted in the jet models, an open mind was kept in re- 
gards to the value of /3. Preliminarily, a relativistic jet was 



2 DESCRIPTION OF THE JET MODEL AND 
THE FIT TO THE DATA 

In this section, the parameter s in the jet model 
are defined. One should review iGhisellini et al (|l985l ): 
iGhisellini and Marasch 1 (|l989 l) for further discussion. The 
jet is parameterized in cylindrical coordinates, (p, (p, z). The 
jets models are axisymmetric and have no transverse gradi- 
ents (this assumption is modified below). The jet shape and 
axial gradients are parameterized by, 
7min and 7max through the the expressions. 



p = az"; 20 < 2 < 2max, A''(7) = K-^ ^, 



(1) 



In equation (1), the electron energy is E — ^m^c? and 
the particle number density is A'' = J N{-f)d'y. The mag- 
netic field strength is B. The jet begins and ends at zq and 
2max, respectively. The models typically have two regions. 
The base of the jet is called the paraboloid, which is the 
stronger invers e Com pton source, and it is characterized in 
IGhisellini et all (| 19851 ) by e = 0.5. The outer part of the jet 
is called the conical region typically with e = 1. This is gen- 
eralized to other values of e for the outer jet in the following. 
The outer jet is the primary source of the radio emission. 
Tables 1 and 2 list the jet model parameters and physical 
quantities, respectively, that characterize a jet solution. The 
dimensional input parameters and power law indices in Ta- 
ble 1 define the model in terms of equation 1. The chosen 
value is either fitted to the data or set by the dimensions of 
the physical system as noted in the table. All the quantities 
in Table 2 are physical characteristics of the jet model. It is 
noted in the table that the values are either an input param- 
eter fitted to the data or a derived physical quantity that 
is computed from the model. The jet plasma is protonic. 
Notice that Tables 1 and 2 list four shells and a core (shell 
0) instead of one jet. The reason for this necessity was the 
surface brightness distribution in the VLBA images is very 
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Figure 2. The model fit to the broadband SED is in the left frame. The IR error bars represent the variability at the time of observation 
from Fuchs et al (2003). The radio data is from a simultaneous VLA measurement at 8.3 GHz and a 15 GHz Ryle telescope observation 
2.5 hours earlier. The error bars represent the variability seen in the radio measurements over a two hour time scale (the length of the 
VLBA observation). This plot includes contributions from the oppositely directed, but otherwise identical counter-jet. The model has 
an excess above a simple power law in the sub- mm as in the epochs observed in Ogley et al (2000). The right frame is a close-up of the 
spectral energy distribution in the X-ray region. The INTEGRAL data is in green and the RXTE data is in black. The blue curve is the 
disk blackbody spectrum plus the EIC steep powerlaw from the core. The red curve also includes the EIC from shells 1 to 4 



extended. A jet with no transverse gradients that accurately 
models the SED always produces a surface brightness distri- 
bution that decays too rapidly with z to represent the VLBA 
maps (see Section 4) . Figure 1 shows the relative scale of the 
nested shells near the base. A continuous transverse gradi- 
ent would provide a smoother solution. However, the simple 
crude model of 5 uniform (there are no transverse gradients 
within each zone) zones is shown to be a reasonably accurate 
representation of the spectrum and the brightness profile in 
the following. 

Each shell is parameterized in the tables. Columns 2, 
3, 4, 10 and 11 of Table 1 define the parameters for the 
base of the jet in terms of equation (1) for each shell. Sim- 
ilarly, columns 5-9 are the outer jet parameters. The 
initial data surface from which the jet solution begins is 
formally defined for each shell "i" as the set of points 
Pbase(i - 1) < p < pbaso(«) such that z = 2o(base), where 
zo(base) is the minimum axial coordinate in shell "i" and 
Pbaso(j — 1) = f or i = 0. One impor tant modification 
to the procedure of iGhisellini et all l|l985l ) is that the inner 
shells must transfer radiation not only through the optical 
depth provided by its own synchrotron self absorption, but 
also that of the surrounding shells. This affects the turnover 
frequency between the optically thick region of the spec- 
trum and the op tically thin region of the spectrum through 
equation (11) of IGhisellini et al (1985(). Because the opac- 
ity dies off with cylindrical radius, this is a second order 
correction to the spectrum from each shell, but is retained 
to improve the accuracy of the calculations. In the follow- 
ing models, = 70° is assumed as was determined to be 
consiste nt with the kinematics of disc rete relativistic ejec- 
tions in iMirabel and Rodriguej l|l994h . The flux densities 
are converted to luminosit y using a distance of 12 . 5 kpc to 
the source as estimated in Mirabel and Rodriguez] (Il994l ). 

The baseline fit to the broadband spectral energy distri- 
bution (SED) is in Figure 2 with the X-ray portion magnified 
in the r ight hand fr a me. T hese plots are the analog of Figure 
5 from iFuchs et all (|2003l ). The X-ray data were calibrated 



and reduced with current software suite, HEAS OFT V 6.10 
and OSA V 9.0, that are improvements from iFuchs et all 
(|2003l . [ 20041 ). The new results were generously provided by 
Jerome Rodriguez. The fit in the X-ray region is not that 
good around 10^^ Hz. There is increased curvature between 
6 X 10" Hz and 10^^ Hz that is more evident in the new re- 
duction of the data than in iFuchs et all (|2003h . There is an 
overlap of the data in the HEXTE and the INTEGRA L spec- 
trum in this region as discussed in lFuchs et all (|200>j). Even 
though the two data sets do not agree in this area (the IN- 
TEGRAL spectrum i s flatter wit h high er flux levels, see the 
related discussion in iFuchs et all (|2003l )). this does not ap- 
pear to be the origin of the spectral curvature, but coinciden- 
tal. This region is not well represented by the EIC powerlaw 
from the jet core (shell 0) in blue and is actually exacerbated 
by the spectral hardening that seems to be real at frequen- 
cies above 10^^ Hz. In the model, the hard excess is fit by 
the red, total luminosity curve which includes a second EIC 
contribution from the surrounding shells. In this decompo- 
sition, the spectral curvature at frequencies just below 10^^ 
Hz is consistent with spectral ageing from Compton losses 
of the higher energy electrons in the cor e jet which have the 



e .let \ 

shortest radiative time scale, r ~ 7""^ (jTucker 1975). This 
suggest s that a more elaborate jet m odel is indicated for the 
core. In lRodriguez et all l|2004l2008l ), the poor powerlaw fit 
near 10^^ Hz is improved by including an inverse Compton 
coronal component with a thermal cutoff. The extra param- 
eters introduced by a corona with cutoff is beyond the scope 
of this paper, which is an exploration of how well a simple 
jet model can explain the data. 

Figure 3 is the surface brightness distribution based on 
the modeled jet emissivity as it would appear with the res- 
olution provided by the synthetic beamwidth of the VLBA 
as modified by the wavelength dependent interstellar scat- 
tering {A9 ~ 0.15(A /lcm)^ mas) that was determined in 
iDhawan etall (|2000l ). The plasma in the outer shells pro- 
duce the low frequency emission from innermost part of the 
jet, "the radio core" (or peak flux per beamwidth) and the 
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observed base of the radio jet (the partiaUy resolved emis- 
sion located within one beammwidth of the core, 33.6 AU 
at 8.3 GHz and 17.0 AU at 15 GHz) that is depicted by 
the brightness distribution in Figure 3. Shell 1 produces the 
farthest reaches of the jet. The lowest contour is at 15 times 
the RMS noise of the radio images, thus the vertical error 
bars are very small on all the data points. The dominant 
source of error is the sparse u-v coverage (the missing 8.3 
GHz flux compared to the VLA measurement noted in the 
figure caption). There is really no way to attach these errors 
to any of the individual contours. All that one can say is 
that the relatively long baselines of the VLBA favors missed 
flux on the least compact scales, the outer-most contours. 



3 PHYSICAL INTERPRETATION OF THE 
JET MODEL 

The model presented above and described by the input pa- 
rameters in Table 1 and Table 2 can be interpreted phys- 
ically in terms of quantities like mass flux, Poynting flux 
and magnetic dominance. In the following, it is described 
how the model parameters relate to physical quantities in 
Table 2. Column 2 of Table 2 gives the energy flux in each 
shell at the base. The next 4 columns provide the number 
density, magnetic fleld morphology and component strength. 
Columns 7 and 8 are the ratio of magnetic energy flux to 
mechanical energy flux at the minimum value of z for each 
jet region determined from equation (1) and equations (2) 
-(5), below. The ninth column is the poloidal magnetic flux 
in the shell determined from the same equations. The last 
three columns define the electron energy power law from 
equation (1) that applies throughout the shell. 

The innermost shell labeled "0" (the core) is initiated 
in the innermost part of the accretion disk. In the model, it 
resides in the ergospheric portion of the disk inside of 1.7 M 
(M = 2 X 10^ cm in g eometrized un i ts for the 14 Mq black 
hole mass inferred bv ICreiner et aj (|200ll 'l). Thus, one re- 
quires a rapidly spinning black hole and a black hole spin of 
a/M > 0.99 is necessary to provide sufficient surface area in 
the disk. If a slowly spinning black hole is used in the model, 
the volume of the base of the jet increases dramatically and 
so does Q (which is highly dependent on the total amount 
of energized particles and magnetic field - large jet volumes 
create large total plasma energy in the models). The mag- 
netic field is turbulent and it is assumed to be advected with 
the jet as in iBlandford and Koenigll (|l979l '). The jet begins 
near equipartition as noted in Table 2. From equation 1, the 
value of mbasc = 2.64 means that magnetic energy is dissi- 
pated rapidly in the plasma near the base of the jet. The 
lost energy is the X-ray emission forming the high energy 
power law. The core is the most powerful portion of the jet 
and most of the energy is lost to X-rays by inverse Comp- 
ton c ooling. The X- r ay spectrum is computed per the meth- 
ods oflTuckeil (|l975l '): lGhisellini et all l| 19961 ) : ISambruna et all 
j 19971 ). The X-ray power law spectral index, (p — l)/2, fixes 
p = 5.6. This is a very steep value that is consistent with 
the assumption of intense Compton cooling. The disk emis- 
sion is fit by a simple blackbody (instead of a parameterized 
accretion disk model) with T — 1.5 x lO'^ ° K and a surface 
area on each disk face of 2.4 x lO'^^cm^ (the surface area of 
a thin Keplerian disk for r ^ 5M and a/M = 0.99). The 



spectral fit has ~ 4.0 x 10"^* ergs/sec in thermal disk lumi- 
nosity and ~ 2.5 x 10"^* ergs/s ec in jet lum i nosity (note the 
total matches that inferred bv lFuchs et aj (|2004l ) when the 
fiux overestimate from RXTE/PCA of 16% that they note is 
taken into account). The large EUV peak in the synchrotron 
spectrum in Figure 2 is from the same electrons that pro- 
duce the EIC. The outer portion of the core jet at large z is 
inertially dominated and is a weak synchrotron source. 

Shells 1 to 4 are Poynting fiux dominated jets with a 
field line angular velocity, Q, f , equal t o the relativistic Keple - 
rian velocity, f^Kcp, at the foot point [liightman et all (|l975h . 
The field is organized and satisfies the perfect MHD con- 
dition similar to the electrodynamic models of iBlandfordI 
(| 19761 ) : lLoveiac3 l| 19761 ). 



pBP{d<j)/dt-nF) 



(2) 



The jets have an artificial simplification that the poloidal 
bulk velocity is constant, = 0.4c. From equations (1) and 
(2), at large z, the toroidal magnetic field, B"^ , is much larger 
than the poloidal magnetic field, B^ , and the particle angu- 
lar velocity, d(f)/dt, is small compared to Of (bounded ve- 
locity and angular momentum). Perfect MHD requires that 
the total poloidal magnetic fiux is conserved in each shell, 



B^dAi 



(3) 



where dA±_ is the cross-sectional area in the shell. Equations 
(1) - (3) combine to give a relationship between B^ at the 
base of the jet and B"^ at the start of the outer jet, 



-'-'base ' 



zq (outer) 
PbascfiKcp / V 20 (base) 



(4) 



This relationship allows one to compute (as in Table 2) $ 
and the poloidal Poynting flux, , in both regions of the 
jet, where 



1 

4-K 



^FpB'^B^dAj 



(5) 



The base of the shell 1 jet produces the bulk of the 
synchrotron IR emission. The jet emission in Ks band was 
estimated to be about 40% of the total luminosit y (the bulk 
being from the donor star) in lFuchs et all l|2003l ) consistent 
with the fit in Figure 2. The value of p — 3.3 for the shells 
is chosen to reproduce the noticeable hardening of the X- 
ray spectrum above 10^^ Hz seen in the INTEGRAL data 
in Figure 2 and discussed above. More than one half of this 
excess over the steep core powerlaw arises from EIC at the 
base of the shell 1 jet. The other shells provide the remainder 
of the hard EIC excess and the inner part of the brightness 
profile of the jet in Figure 3. Shell 4 is very energetic because 
the volume is huge - it has a very wide base. Shells 2 and 
3 are required to interpolate the brightness distribution in 
Figure 3 between the "radio core" and base of the observed 
radio jet produced by the plasma in shell 4 and the distant 
regions of the brightness distribution that is emitted from 
plasma in shell 1. 

In summary, if one assumes that the hard X-ray power 
law is produced in the jet (as in the fiducial model described 
in the last two sections) then the results in Table 2 imply 
that the jet power is Q = 4.2 x 10^® ergs/sec. Assuming 
a/M = 0.99, the mass fiux in the jet is 7.8 x 10^'^ g/s. 
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Figure 3. The brightness profiles at 8.3 GHz and 15 GHz obtained from the contours of the VLBA images in Fuchs et al (2004). Each 
dot represents a contour level. For comparison, the synthesized flux per beam that is produced from the 5 shell jet model convolved 
with the VLBA beamwidth and interstellar scattering is plotted. The diffuse nature of the jet emissivity requires a very broad base and 
transverse gradients in the jet, greatly increasing the volume and the required power. The excess flux at large z provided by the model 
is a desired effect, since simultaneous VLA measurements at 8.3 GHz in Fuchs et al (2003) show that > 20% of the flux is resolved out 
by VLBA and likely resides at large distances, > 10^^ cm. The plots do not cover the counter-jet as it was noted in Ribo et al (2004) 
and in the text that a single velocity cannot produce the asymmetry in the velocity profile, so there cannot be a good fit to this data 
simultaneously. This might be evidence that transverse and axial gradients in the bulk velocity are required or the jet model is not 
representative of the emission. 



Table 1. Input Parameters of the Jet Component Models 



Shell rribaso "base "base "loutcr riouter doutor e 2o(outer) ZQ (base) Phase 

cm cm cm 

aba ab aaa a c 






2.64 


2 


1.52 X 10^ 


2.4 


2 


3.82 X 10^ 


0.45 


1.0 X 10* 


5.0 X 10*^ 


3.4 X lO" 


1 


1.08 


2 


4.02 X 10^ 


2.1 


2 


2.20 X 10^ 


0.595 


1.5 X 10^^ 


5.0 X lO** 


8.8 X 10^ 


2 


1.02 


2 


6.03 X 10^ 


1 


2 


1 


0.8 


4.0 X 10^2 


5.0 X 10*^ 


1.4 X lO'' 


3 


0.99 


2 


1.35 X 10* 


1 


2 


5.0 X 10-2 


0.92 


8.5 X 10^2 


5.0 X lO'' 


3.0 X lO'' 


4 


1.00 


2 


3.57 X 10-* 


1 


2 


1.4 X 10-2 


1 


6.5 X 10^2 


5.0 X 10^ 


8.0 X lO'' 



Table 2. Table notes: a. input parameter for model fitted to SED data, b. set by baryon conservation, c. input parameter is a minimum 
value consistent with the black hole dimensions 



4 PARAMETRIC SENSITIVITY 

It is not claimed that this fiducial model is the unique jet so- 
lution that fits the data. This raises the question of how the 
model depends on the parametric assumptions. Although a 
formal parametric sensitivity analysis is beyond the scope of 
this work, many trial models were produced and evaluated 
in the process of arriving at the model presented in the pre- 
vious sections. Here are some of the observations that were 
made. 

• Number of Shells The integrated kinematic quanti- 
ties such as the magnetic, energy and mass fluxes are not sig- 
nificantly affected by the number of shells. The more shells, 
the smoother the brightness distribution and the better the 
fit is in Figure 3. For example, if 4 or 6 shells were used in- 
stead of 5, the SED can be reproduced. However, the fit in 
Figure 3, would be considerably worse for 4 shells and some- 
what better for 6 shells. The luminosity in each shell decays 
on scales that are comparable to the distance, z, that the 
peak of the luminosity is from the accretion disk. For exam- 
ple, the 15 GHz luminosity is produced at smaller z than the 
8.3 GHz luminosity in each shell. Thus, the surface bright- 



ness dies off with a shorter scale length at 15 GHz than at 8.3 
GHz. Compare the right hand frame of Figure 3, to the left 
hand frame - the 15 GHz fit is worse in the sense that there 
are large curved sections between actual data points and a 
pronounced infiection point neither of which exist within the 
sparsely samp led data (there are only 6 contour levels in 15 
GHz image in iFuchs et all l|2004l )). This is manifested as the 
almost exponential decay near 10 AU, followed by a sharp 
peak at ~ 25 AU. The more shells, the less pronounced the 
alternating regions of exponential dips and sharp peaks be- 
come. With 4 shells, the 8.3 GHz plot would have similar 
large deviations from the brightness profile as are seen in 
the existing 15 GHz plot in Figure 3. 

• Energy Flux Estimate Given the assumptions in the 
model (in order of importance: the hard X-ray luminosity is 
from the jet, the hard X-ray excess above 10^^ GHz is real 
and from the jet, and the IR excess above the stellar contri- 
bution is from the jet), the model is very close to a minimum 
energy flux. This is obvious since the energy flux is only 1.5 
times the radiated luminosity and the jet power must ex- 
ceed the broadband luminosity. The energy flux grows very 
rapidly, if the jet is assumed to originate at large cylindrical 
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Table 3. Physical Parameters of the Jet Component Models 



Shell 


Q 


N 


B Field 




BP 


5/KE 


5/KE 


$ 


P 


^min 


'7m ax 




10'^* ergs/sec 






10^ G 


106 G 


Base 


Outer 


G - cm^ 










b 


a 


a 


b 


b 


b 


b 


b 


a 


a 


a 





3.37 


9.5 X 10^^ 


Turbulent 


29.8 


29.8 


2.3 


1.7 X 10^2 


Net = 


5.6 


1.5 


1000 


1 


0.12 


2.9 X 10" 


MHD Jet 


6.9 


3.3 


4.8 X 10'' 


2.1 X 10"* 


1.4 X 10^1 


3.3 


50 


1000 


2 


0.04 


1.3 X lO^o 


MHD Jet 


2.6 


2.3 


3.4 X 10^ 


2.6 X 10^ 


8.4 X 10^" 


3.3 


50 


1000 


3 


0.18 


2.7 X 10^ 


MHD Jet 


2.1 


3.0 


1.1 X lO" 


1.1 X lO** 


4.7 X 10^1 


3.3 


50 


1000 


4 


0.49 


8.4 X 10* 


MHD Jet 


1.1 


2.9 


1.0 X 10^ 


1.0 X lO'' 


2.0 X 10^^ 


3.3 


50 


1000 



Table 4. Table notes: a. input parameter for model fitted to SED data, b. derived quantity from parameters of the jet model 



radius. The volume of plasma and fields increases rapidly 
and the energy flux quickly rises above 10"^^ ergs/sec as the 
cylindrical radii of the base of the shells increase by a factor 
of >2. 

• Number Density The number of energetic particles 
is fixed by the location of the base of the jet, the thermal 
X-ray flux from the accretion disk and the constraints on 
the observed external Compton luminosity. The particle flux 
and energy flux are concentrated at smaU radius in the core. 
Driving most of the energy and mass flux from the inner re- 
gion of the disk is favored on physical grounds. It co-locates 
the most energetic shell and largest mass flux at the largest 
potential source of power to drive the jet - the most ener- 
getic part of the accretion flow (fastest rotation rate and 
highest temperature). Moving the core to the smallest pos- 
sible cylindrical radius also moves the other shells inward 
to smaller radius and larger Keplerian velocity, thereby de- 
creasing the total volume and magnetic flux required to drive 
the flow, by equation (5). This minimizes the energy flux in 
the outer shells. In this decomposition, the mass flux in the 
outer shells is constrained by the magnitude of the hard X- 
ray excess above lO'^® Hz (EIC from shells 1 -4) discussed 
in section 2. The large radio fluxes from the outer shells 
combined with a modest EIC X-ray component result in the 
magnetic dominance for the outer shells depicted in columns 
7 and 8 of Table 2. This analysis does not preclude the pos- 
sibility of a different transverse distribution of mass flux in 
a jet that can reproduce the data, however the energy flux 
is likely to be larger. 

• Magnetic Field Strength The magnetic fleld 
strength, B is set by the size of the base of the jet, the 
particle number density and the synchrotron portion of the 
SED. The magnetic field must also be of sufficient strength 
to be the ultimate energy reservoir that supports the bulk 
of radiation losses as discussed in section 3. 

• Electron Energy Spectrum The electron spectrum 
is fixed by the X-ray power law spectrum as mentioned in 
section 2. The only parameter that is not affected by the 
data is Vmax- Since the energy distributions are steep, there 
are very few high energy electrons. A value of Tmax ~ 1000 
was chosen arbitrarily. Any value larger than 700 will show 
negligible differences in the results. 

• Axial Dimensions The value of zo{hase) is set by the 
physical dimensions of the accretion/disk black hole system. 
The quantities defined at zo(base) form the initial data for 
the initial surface of the entire jet. By contrast, 2o (outer) is 
set by the SED fit. For example, large values of zo (outer) 
over-produce the high frequency synchrotron emission. By 



definition (outer) = 2niax(base). The solution is insensitive 
to 2:max (outer) = 5 X 10^^ cm, since the jet emissivity is very 
low in this region. 



5 ALTERNATIVE MODELS 

The shell decomposition in Table 2 is very useful for explor- 
ing variants of the jet model. The most obvious variant is 
to assume that the disk corona and not the jet is the source 
of the X-ray powerlaw. By lowering the mass density and 
magnetic field density significantly (say, a factor of 100 for 
the sake of argument) in the core, one can simply ignore the 
energy flux and luminosity contribution to the composite 
jet. Then summing the remaining components in Table 2, 

Q ~ 8.3 X lO^^ergs/sec , for a coronal X — ray source, (6) 

and the jet plasma still reproduces the SED (less the coronal 
X-rays) and the radio brightness profile in Figures 2 and 3, 
respectively. 

There is uncertainty in the non-relativistic jet velocity 
as noted in Section 1. Thus, it desirable to know the sensitiv- 
ity of the jet power to different values of jet poloidal velocity, 
Q{vP). The shell model lends itself to estimating Q for dif- 
ferent values of . First decompose Q into 3 discrete pieces. 
Then note from from equations (2) and (5), that in the 
outer jet has a simple scaling with , Souter ~ [B'^'Y . 
Since Q is a sum of radiation losses, the mechanical kinetic 
energy flux in the outer jet, KE, and Soutor, one can use 
the results in Table 2 to parameterize Q as a function of 
when the X-ray power law is from t he jet. Implementing t he 
Doppler enhancement factors from iLightman et all l|l975h . 



S{I3 = 0.4, e = 70°f + 6{I3 : 



-0.4, 9 = 70° 



<S(/3, e = 70°)* + 6{-P, 9 = 70°)* 



x(2.5 X 10^^) 



+ 



-0.4, 9 = 0°) 



(F - l)(«^/0.436c)(8.5 x 10^^) 



-H(«^/0.4c)(8.3 X 10^^)ergs/sec, (7) 

where S = S{I3, 9), the Doppler factor that was introduced 
in Section 1, is a two parameter function of /3 and 9. In this 
notation, (5(/3 = 0.4, 9 = 70°) means to evaluate the Doppler 
factor at /3 = 0.4 and 9 = 70°. The jet poloidal velocity is 
yP = cp. 

The Doppler correction on the flrst term in equation (7), 
means that larger 5 factors require smaller intrinsic radia- 
tion losses to produce the observed X-ray powerlaw SED in 
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Figure 2. The Doppler correction on the second term resuhs 
from the redshifting of the disk thermal emission in the rest 
frame of the plasma at the base of the jet. A slow jet sees a 
stronger seed field and requires fewer hot particles to repro- 
duce the EIC SED and therefore has a smaller bulk kinetic 
energy flux. The last term is the electromagnetic energy flux 
in the outer jet. As noted in section 1, this formula will not 
hold for large j3 since SSC contributions become important. 

Similarly, if the disk corona is the source of the X-ray 
powerlaw, one can use the results in Table 2 and equation 
(6) to parameterize the jet power as a function of jet velocity 
as 



Q{v^) ^ [u^/0.4c][8.3 X 10^'^] ergs/sec 



6 DISCUSSION 



(8) 



The crude model presented here has demonstrated a few 
facts about the compact jet in GRS 1915-1-105 in the "high 
plateau state": 

• It is possible to model the jet in a manner analogous to 
AGN 

• The jet is powered during a high accretion state with a 
thermal luminosity of ~ 4 x 10^* ergs/sec, corresponding to 
an Eddington rate of « 23%. 

• The jet brightness profile requires a highly stratified jet 
in the transverse direction 

• Relativistic jet bulk velocity is not favored based on the 
huge implied energy flux, > 3.5 x 10*'' ergs/sec if the jet is 
the source of the X-ray power law. 

An interesting consequence of the model was that 
a rapidly spinning black hole was preferred on energetic 
grounds. This was not to tap black hole spin energy, but 
to use the surface area and high angular velocity of the 
plasma in the ergosphere to power the jet core. It is interest- 
ing to note that the 3- D, perfect MHD, general r elativistic 
numerical simulations in lHawlev and KrolikI l|2006l ) find very 
high efficiency "funnel wall jets" emanating from this region 
of the ergosphere of high spin black holes. The outfiow is 
somewhat similar to what is described in the core shell in 
this treatment, near equipartition with turbulent disordered 
magnetic fiux and an outflow velocity of ~ 0.3c. The effi- 
ciency is sufficient to power the core if a/M > 0.95. They ar- 
gue that the power source is pressure gradients driven by the 
disk which seems energetically reasonable in GRS 1915-1-105 
considering the large thermal luminosity noted above. 

Even though there is some numerical evidence to sup- 
port the notion of a near-relativistic jet core (where most of 
the energy transport occurs), it is not clear if there is any 
analog of the compact jet in the AGN population. In par- 
ticular, are there near-relativistic outflows that are highly 
stratified in the transverse direction that carry a mass flux 
comparable to the accretion fiux? It is also not obvious if 
there is any structure that can be construed as an analog of 
a non-relativistic high luminosity, optically thick, radio jet 
in AGN. If there is scale invariance in black hole accretion 
systems then there should be an analog of the compact jet 
in some AGN. The only obvious way to reconcile this dif- 
ference is to consider the possibility that an AGN radio jet 
analogy is not appropriate and the analog of the compact jet 
in AGN takes a different guise due to the lower accretion disk 



temperature. The AGN feature that possesses the majority 
of the characteristics of the jet model of GRS 1915-1-105 
are the near-relativistic outfiows that are seen in X-ray ab- 
sorption. A few radio quiet quasars have been observed to 
have out-flowing. X-ray a bsorbing winds with velocities be- 
tween 0.1 c and 0.6 c (IChartas et all [20031 : iGofford et all 
I2OIII : iPounds et al l2003l : ISaez et all 120091 ') ~The outflow in 
APM08279-I-5255 is the most energetic with a mass flux 
similar to the accretion mass flux and a kinetic luminosity 
comparable to the the rmal luminosity of the accretion flow 
as in GRS 1915-1-105 l|Saez et al|[2009l ). The fastest known 
out-flowin g, X-ray absorbing w inds in radio loud AGN have 
V < 0.15c (|Tombesi et alll2010^ . 

In summary, due to the relatively low jet velocity and 
the large transverse stratification, perhaps the compact jet 
in black hole binaries is best described as a bipolar near- 
relativistic wind. In this analogy, it is not obvious why the 
Galactic black hole version of the wind is luminous in the 
radio band and the AGN version is not. 

It has been argued in iKoerding et all (|2006l ) that there 
may exist a compact, optically thick jet in low- luminosity 
AGN. This may or may not be applicable to the study here 
in that the thermal luminosity was 23% of the Eddingtom 
limit in these observations of GRS 1915-1-105. In any event, 
it is worthwhile to study jets from low luminosity AGN with 
VLSI to look for evidence of nonrelativistic optically thick 
winds or jets. 

The implication of these jet models to other Galac- 
tic black holes can only be speculated upon because other 
sources do not have the wealth of high resolution radio im- 
ages. Furthermore, GRS 1915-1-105 is a very strong X-ray 
source, so there are many excellent studie s of the radio 
j et/ X- r ay connection. Consider , the results of lRushton et all 
l|2010bD : lRodriguez et all (|2008 h that the existence of a hard 
power-law component is tantamount to the existence of a jet. 
Within the framework of the assumption that inverse Comp- 
ton emission in the jet is the dominant source of the power 
law component, the jet power can be estimated from the 
model presented here as approximately 1.5 times the hard 
X-ray power law luminosity. This estimate can be applied 
to other Galactic black holes in the hard state. However, it 
would be much more accurate to constrain the energy flux 
derived from the basic jet model on a case by case basis 
with deep radio images and broadband simultaneous spec- 
troscopy, which are not readily available for other Galactic 
black holes. There are estimates of the compact jet energy 
flux based on less direct methods for other black hole candi- 
dates such as Cygnus X-1. For example, the emission from 
a nebular ring 5 pc from Cygnus X-1 has been used to esti- 
mate the energy flux required to energize the structure. This 
in turn is related to the compact jet power on much smaller 
scales by assuming an invisible flow of energy to large scales. 
This leads to estimates of th e kinetic power o f the compact 
jet o f lO^'' to lO^'^ ergs/sec (|GaUo et alll2005l : iRusseU et all 
[2003). This estimate is 0.3 to 1.0 of the total radiative lu- 
minosity, which is similar to the jet power of approximately 
0.6 times of the total radiative luminosity that was found 
for the compact jet in GRS 1915-1-105 here. 

Another consequence of the model is that a relativistic 
outflow velocity in compact jets is not favored based on ener- 
getic concerns. It is interesting to compare this to less direct 
estimates of the compact jet speed in GRS 1915+105 and 
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other Galact ic black hole candidates . Usin g the very basic 
jet model of iBlandford and Koenigll ([l9zi), it is argued in 
ICasella et all (|2010l ) that r > 2 in GX 339-4, which was have 
ruled out for GRS 1915+105 in the epoc h considered here. A 
very g eneral argument was presented in lHeinz and Merlonil 
l|2004 ) that was based on total X-ray and radio fluxes from 
a one zone model (that happens to be cylindrical in shape). 
They found a loose set of constraints that favored F > 2 
in Galactic X-ray binaries, but w as also con s istent with 
j3 ~ 0.5. Contrary to these works iGallo et" 

3 (l2003l ) used 

Monte Carlo simulations to argue that F ^ 2 is preferred for 
Galactic X-ray binaries in the low/hard state, marginally 
consistent with the energetic constraints of the models pre- 
sented here. 
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